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The next to the leading order Casimir effect for a real scalar field, within 4> 4 theory, confined 
between two parallel plates is calculated in one spatial dimension. Here we use the Green's function 
with the Dirichlet boundary condition on both walls. In this paper we introduce a systematic 
perturbation expansion in which the counterterms automatically turn out to be consistent with 
the boundary conditions. This will inevitably lead to nontrivial position dependence for physical 
quantities, as a manifestation of the breaking of the translational invariance. This is in contrast to the 
usual usage of the counterterms, in problems with nontrivial boundary conditions, which are either 
completely derived from the free cases or at most supplemented with the addition of counterterms 
only at the boundaries. We obtain finite results for the massive and massless cases, in sharp 
contrast to some of the other reported results. Secondly, and probably less importantly, we use a 
supplementary renormalization procedure in addition to the usual regularization and renormalization 
programs, which makes the usage of any analytic continuation techniques unnecessary. 

PACS numbers: 



I. INTRODUCTION 

During the last fifty years many papers have been writ- 
ten on the calculation of the Casimir energy. In this 
paper we introduce a new approach in regards to the 
renormalization program. We found it most suitable to 
introduce our approach in the simplest example possi- 
ble, i.e. a real scalar field confined between two parallel 
plates in 1+1 dimensions, with <j> self- interaction. As we 
shall see, our results for the next to leading order term 
(NLO) are finite for both the massive and massless cases 
and this differs significantly from what exists in the lit- 
erature. It is therefore suitable to start at the beginning. 
In 1948 H.B.G. Casimir found a simple yet profound ex- 
planation for the retarded van der Waals interaction [l| . 
After a short time, he and D. Polder related this effect to 
the change in the zero point energy of the quantum fields 
due to the presence of nontrivial boundary conditions [2| . 
This energy has since been called the Casimir energy. 
The zero-order energy in perturbation theory has been 
calculated for various fields (see for example 3). Also 
the NLO correction, which is usually called the first-order 
effect, has been computed for various fields. For the elec- 
tromagnetic field this correction is said to be due to the 
following Feynmann diagram ^Q^) , and has been com- 
puted first by Bordag and collaborators 
However, note that this correction is a two loop correc- 
tion in this case and is C(e 2 ). Moreover the two-loop ra- 
diative corrections for some effective field theories have 
been investigated in [l(l EH, 03- Next, in the case of 



a real massive scalar field NLO correction to the ener. 
has been computed in 



\A 14, 15, \m 17, 18, 19, \2 



This correction is a two loop correction in this case but 
is 0(A). Moreover, N. Graham et al used new approaches 
to this problem by utilizing the phase shift of the scat- 
tering states 21|, or replacing the boundary conditions 



by an appropriate potential term [22|. However, the au- 
thors use the free counterterms, by which we mean the 
ones that are relevant to the free cases with no nontrivial 
boundary conditions, and are obviously position indepen- 



dent. Only in Ref. 17| the author notes that in certain 



cases, counterterms can depend on the distance between 
the plates. The first use of nontrivial boundary condi- 
tions for the renormalization programs in problems of 
this sort seems to be due to Fosco and Svaiter [23[ . These 
authors use free counterterms in the space between the 
plates and place additional surface counterterms at the 
boundaries. Later on various authors proposed the use of 
exactly the same renormalization procedure for various 
physical problems The first calculation for the NLO 
of Casimir energy for the massive scalar field using this 
renormalization program is done in Ref. 25J . We should 



note that their results for the massless limit in 1+1 di- 
like those of 



mcnsions, 



19j, is infinite. It is also worth 



mentioning that all the papers on the analogous calcula- 
tions of the NLO corrections to the mass of solitons, that 
we are aware, of use free counterterms (see for example 



28(). In references [27J the authors used the mode 
number cutoff introduced by R.F. Dashen (1974) 28j to 
calculate the NLO Casimir energy due to the presence of 
solitons. 
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In this paper, we present a systematic approach to 
the renormalization program for problems which are 
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amenable to renormalized perturbation theory, and con- 
tain cither nontrivial boundary conditions or nontrivial 
(position dependent) backgrounds, e.g. solitons, or both. 
Obviously all the n-point functions of the theory will have 
in general nontrivial position dependence in the coordi- 
nate representation. This is one of the manifestations of 
the breaking of the translational symmetry. The proce- 
dure to deduce the counterterms from the n-point func- 
tions in a renormalized perturbation theory is standard 
and has been available for over half a century. Using 
this, as we shall show, we will inevitably obtain posi- 
tion dependent counterterms. Therefore, the radiative 
corrections to all the input parameters of the theory, in- 
cluding the mass, will be in general position dependent. 
Therefore, we believe the information about the non- 
trivial boundary conditions or position dependent back- 
grounds are carried by the full set of n-point functions, 
the resulting counterterms, and the renormalized param- 
eters of the theory. Our preliminary investigations have 
revealed that our position dependent counterterms ap- 
proach the free ones when the distance between the plates 
is large. Their main difference is for positions which are 
about a Compton wavelength away from the walls, al- 
though it is also nontrivial at other places. Here we use 
this procedure to compute the first-order radiation cor- 
rection to the Casimir energy for a real scalar field in 
1+1 dimensions. We compute this correction for both a 
massive and a massless scalar fields and show that the 
massless limit of the massive case exactly corresponds to 
the massless case. 

In addition, up to now all the papers on the Casimir ef- 
fect, that we are aware of, use some from of analytic con- 
tinuation. We share the point of view with some authors 



such as the ones in [14j, [18| that the analytic continuation 
techniques are not always completely justified physically. 
Moreover, like the first of the aforementioned authors, 
we have found counterexamples, which we point out in 
this paper and elsewhere 29( . The counterexamples show 
that it alone might not yield correct physical results, and 
sometimes even gives infinite results [30(. Therefore, we 
prefer to use a completely physical approach by enclosing 
the whole system in a box of linear size L, which even- 
tually can go to infinity, and calculating the difference 
between the zero point energies of two different configu- 
rations. The main idea of this method is actually due to 
T.H. Boyer 3l|, who used spheres instead of boxes. This 
we shall call the "box renormalization scheme" and can 
be used as a supplementary part of other usual regular- 
ization or renormalization programs. This box renormal- 
ization scheme, has the following advantages: 

1. Use of this procedure removes all of the ambigui- 
ties associated with the appearance of the infinities, 



and we use the usual prescription for removing the 
infinities in the regulated theory, as explained in 
Sec. MI Bl This is all done without resorting to 
any analytic continuation schemes. 

2. The infrared divergences which generically appear 
in these problems in 1+1 dimensions automatically 
cancel each other. 

3. In order to calculate the Casimir energy we sub- 
tract two physical configuration of similar nature, 
e.g. both confined within finite regions, and not one 
confined and the other in an unbounded region. 

4. This method can be used as a check for the cases 
where analytic continuation yields finite results, 
and more importantly, can be used to obtain finite 
results when the former yields infinite results. 

we should mention that some authors believe that use 
of box regularization or renormalization procedures, in 
which the size of the box eventually goes to infinity could 
be avoided by using appropriate boundary conditions on 
the fields at spatial infinity [32j 

In Section [TT| we calculate the leading order term for 
the Casimir energy in d space dimensional case. We do 
this first of all to explain more completely the physical 
content of the problem and set up our notations. Sec- 
ondly this computation is just about as easy as to do in 
d dimensions as is in the one dimensional case. In Section 
Mil we compute the first order radiative correction to this 
energy. In order to do this we first state the renormaliza- 
tion conditions, and then derive expressions for the the 
first order radiative corrections for both the massive and 
massless cases. We show that the results for the mass- 
less limit of the massive case and the massless case are 
equal. In Section [TV1 we summarize our results and state 
our conclusions. 



II. THE LEADING TERM OF THE CASIMIR 
EFFECT 

The lagrangian density for a real scalar field with <f> 
self-interaction is: 



1 



C(x) = -[d^{x)Y 



1 2 I \2 ^0 , x4 

-m Q tp(x) - -^v{x) , 



(1) 



where m and Ao are the bare mass and bare coupling 
constant, respectively. Here we calculate the leading 
term for the Casimir energy in d spatial dimensions. Ob- 
viously the leading term, in contrast to the higher or- 
der corrections, is independent of the form of the self- 
interaction. The Casimir energy is in general equivalent 
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to the work done on the system for bringing two paral- 
lel plates from ±00 to ±a/2. As mentioned before, part 
of our renormalization procedure is to enclos the whole 
system in a d dimensional cubical box of sides L. To 
compute this leading term, we first compare the energies 
in two different configurations: when the plates are at 
±a/2 as compare to ±6/2. We name the axis perpen- 
dicular to the plates the z axis. To keep the expressions 
symmetrical, we choose the coordinates so that the edges 
of the confining box are at ±L/2 in any direction. 



where, 
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(5) 



Here and k ai _ n denote the momenta parallel and per- 
pendicular to plates (in z-direction), respectively. Also 
aj, and creation and annihilation operators obey- 

ing the usual commutation relations: 

[ a 7i) a n'] — fin,n' j [ a n, a n'] = [ a n> a rc'] = 0; 

and a|0) =0 defines the vacuum state in the presence of 
boundary conditions. Using the above equations one can 
easily obtain 



L 


a 
2 


a 
2 


I 

2 






a2 


a1 


a2 





FIG. 1: The geometry of the two different configurations 
whose energies are to be compared. The labels al, etc. de- 
note the appropriate sections in each configuration separated 
by the plates. 



The total zero point energy of the upper configuration 
in figure ([I]) will be called Eb and of the lower one E a . 
In our box renormalization scheme we need to define the 
Casimir energy as follows 



L d ~\ / dkk d - 2 Yuj ain , 

2 (27r)d-iy ^ 



(6) 



where 7v ' denotes the usual free Hamiltonian density, 
easily obtained from the Lagrangian density, and the su- 
perscript (0) denotes the zero (or leading) order term of 

27T d / 2 

this energy. Also k = Ik^l, and fid = _ . is the solid 



angle in d-dimensions. Therefore, 
L^ 1 fid_i 



r(l) 



E&) _ E (o) = 



2 (2tt 



id— 1 



dkk 



d-2 



>("), (7) 



lim 

6/ a — >oo 



lim (E a -E b ) 

L /b — >oo 



where, 



£-a — E ai + 2E a2 , Eb — Eb 1 + 2Eb 2 



(2) 



(3) 



Here we choose the Dirichlet boundary condition on the 
plates. Then we can expand the field operator ip in the 
eigenstate basis appropriate to this boundary condition, 
and its explicit second quantized form, for example in 
region a\ becomes 



<P ai 0) 




sm 



k ai ,n(z + ) 



k ai ,n{z + — ) 



(4) 



where, 



g(n) = aj ai ,n + 2uj a 



2uj b2 



Now we are allowed to use the Abel-Plana summation 
formula, since we now expect the summand to satisfy 
the strict conditions 40j for the validity of this formula. 
That is, we expect any reasonable renormalization pro- 
gram for calculating any measurable physical quantity to 
yield finite results. The Abel-Plana summation formula 
gives 



E (o)_ E (o) = Ld 1 "rf-i 
6 2 (2iT) d - 1 



-.9(0) 



g(x)dx + i 
o Jo 



dkk d ~ 2 

50 g(it) - g(-it) 



„2irt 



1 



dt. 



,(8) 



where g(0) vanishes in this case due to our box renormal- 
ization. The second term in the bracket, using suitable 
changes of variables, becomes 
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- [°° dn (m 2 + k 2 + k 2 ) 1/2 + 2^— ^ dn (m 2 + k 2 + k 2 ) 1/2 

71 JO 2-7T J 

°l d K (m 2 + k 2 + k 2 ) 1/2 -2^- / dn(m 2 + k 2 + K 2 ) 1/2 = 0, (9) 



Jo ' 2tt 7 

72 7T 

where k for example in the first term denotes — , obviously treated as a continuous variable. The above calculation 

a 

shows that this term is exactly zero. Therefore, only the branch-cut term (the last term in Eq. ([8])) gives nonzero 
contribution and the final result is 

2L d -i m ( n d+1)/2 ^ 1 



p(0) _ F (o) _ ^ " l \ " 

( 47r )(d+l)/2 A^j{d+l)/2 

\ K (.d+iy2( 2a 3 m o) K (d+1)/2 (2bjm ) | 2K {d+1)/2 [(L - a)jm ] 2K (d+1)/2 [(L - b)jm Q ] \ 



a (d-l)/2 5(d-l)/2 1 ^ L- a yd-i)/2 ^ L-b y<i-i)/2 



where K n (x) denotes the modified Bessel function of or- 
der n. Using Eq. @ for the Casimir energy and noting 
that K n (x) is strongly damped as x goes to infinity, only 
the first term remains when the limits are taken, and the 
result is 



E (o) 



2L 



d-1 



(d+l)/2 
t 



(4 7r )(d+l)/2 a (d-l)/2 

.7 = 1 



K 



(d + l)/2 



(2ajm ) 



j(d+l)/2 



If we set d = 3, we have 
E (o) 



L 2 m 2 y, K 2 (2ajm, ) 



8ir 2 a 



3 = 1 



r 



with the following limits, 



E (0) 

Cas. 



-L 2 



1 



-L 2 7T 



8tt 2 o ^ 2a 2 j 4 1440a 3 



as mo 



(11) 
(12) 





-L 2 mp 



3/2 — 2aniQ 



k 8\/2 to 



as amp ^ 1. 



(13) 

The results are in agreement with what exists in litera- 
ture (see for instance 0, EH)- H 1S interesting to note 
that for the massless case, the result is, not surprisingly, 
exactly half of the corresponding expression for the elec- 
tromagnetic case. 

For the d = 1 case, Eq. (fTTj) becomes 



E cL=-^\ " : 



(14) 



and its limits are 



4°2 



7T 

24^ 



as 777-0 



4iPe- 2am » as am »l, 
4 V ira 



(15) 



It is interesting to note that if we solve the massless case 
exactly the branch-cut terms simplify to give 



£(0) 



t 



-dt 



1 



1 

4a7r 



C(2) 



7T 

24rl 



which is identical with our result for the massless limit, 
and is reported for example, in Refs. [jjl ]. 



III. FIRST-ORDER RADIATIVE CORRECTION 

Now we calculate the next to the leading order (two 
loop quantum correction) shift of the Casimir energy for 
a scalar field in 4 theory using the renormalized pertur- 
bation theory in 1 + 1 dimensions. As mentioned before, 
the main idea of our work is that when a systematic treat- 
ment of the renormalization program is done, the coun- 
terterms needed to retain the renormalization conditions, 
automatically turn out to be position dependent. This, 
as we shall see, will have profound consequences. How- 
ever, our main scheme of canceling the divergences using 
counterterms and a few input experimental parameters, 
is in complete conformity with the standard renormal- 
ization approach. To set the stage for the calculations, 
we shall very briefly state the renormalization procedure 
and conditions. 



A. Renormalization Conditions 

The <t> A Lagrangian Eq.JT]), after rescaling the field 
ip = Z 1 l 2 Lp Tl where Z is called the field strength renor- 
malization, and the standard procedure for setting up 
the renormalized perturbation theory, becomes (see for 
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example [34j|). 



= \[d^r{x)] 2 - ^m 2 ip r (x) 2 - ^ r (a;) 4 

(16) 

where S m , 6\, Sz are the counterterms, and m and A are 
the physical mass and physical coupling constant, re- 
spectively. In this problem we are to impose boundary 
conditions on the field at the walls. An alternative ap- 
proach would be to add appropriate external potentials 
to the Lagrangian so as to maintain the boundary con- 
ditions on the fields. We will use the first approach. Ob- 
viously the presence of nontrivial boundary conditions 
breaks the translational invariance and hence momenta 
will no longer be good quantum numbers. Therefore we 
find it easier to impose the renormalization conditions 
in the configuration space. For example, the standard 
expression for the two-point function is, 



(si\T{<i>(xi)4>{x2)}\n) 



lim 

T— >oo(l— ie) 



/ V^{x 1 )^{x 2 )e 



i /5 T cd A x 



}|0) 



JV<P' 



|0) 



(17) 



Since the birth of quantum field theory, as far as we 
know, the assertion has always been that the above ex- 
pressions can be expanded systematically when the prob- 
lem is amenable to perturbation theory. For example, in 
the context of renomalized perturbation theory, as indi- 
cated in Eq. (|16|) . we can symbolically represent the first 
few terms of the perturbation expansion of Eq. (|17p by 



Xi \^ X? 



12 



x 2 xi x x 2 Xl 



-e- 



*2 



(18) 

© * refers to the appropriate counterterm. 



where 

x x x x 2 

It is obvious that the above expression represents a sys- 
tematic perturbation expansion, and most importantly, 
all of the propagators on the right hand side should be 
the one appropriate to the problem under consideration, 
that is they should have the same overall functional form 
as the first term. Our first renormalization condition is 
that the renormalized mass m should be the pole of the 
propagator represented by the first term in (jT5J). This im- 
plies the second and third diagrams should cancel each 
other out in the lowest order, and this in turn implies the 
cancelation of the UV divergences in that order, and that 
the counterterms will in general turn out to be position 
dependent. The renormalized mass m will naturally turn 
out to be position dependent as well. However, we only 



need to fix the value of m(x) at one position between 
the plates by our renormalization condition. The exact 
functional dependence of m(x) will then be completely 
determined by the theory. That is, we insist the overall 
structure of the renormalization conditions such as above, 
and the counterterms appearing in them should be deter- 
mined solely from within the theory, and not for example 
be imported from the free case. The equations are self 
deterministic and there is no need to take such actions. 
Obviously we still need a few experimental input param- 
eters for the complete renormalization program, such as 
m(x) for some x. Analogous expression and reasonings 
could be easily stated for the four-point function. 

To one-loop order the renormalization conditions de- 
rived from Eq. (|18|) and its four-point counterpart, are 



S z (x) = 0, S m (x) 
and S\(x) = 0, 



(19) 



respectively. Here G(x,x') is the propagator of the real 
scalar field and x = (t,z). Obviously the counterterms 
automatically incorporate the boundary conditions and 
are position dependent, due to the dependence of the 
two and four-point functions on such quantities. Now, 
the higher order contributions to the vacuum energy in 
the interval al (i.e. z £ [~^, §]) is 



AE ai = E$ + E£ 



a/2 



-a/2 



dz(n|w z |fi) 




,(20) 



where^g^= — iS\(x) and — — = i[p 2 5z{x) — 5 m {x)) 

refer to the counterterms. Accordingly, the 0(A) contri- 
bution to the vacuum energy is 




~^ m (x)G ai (x, x) 



(21) 



where G ai (x, x') is the propagator of the real scalar field 
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in region al. Using Eqs. (fl9|) and (I2ip . we obtain 

-A r /2 



ai 



-a/2 



G ai (x 7 x)dz. 



(22) 



B. The Massive Case 

As mentioned before, here we choose the Dirichlet 
boundary condition on the plates. Then, after the usual 
wick rotation, the the expression for the Green's function 
in the two dimensional Euclidean space becomes 

2 f dw ujt'-t) \ " Sm [ k ai,n( Z + f )] Sin [ k ai,n( Z ' + §)] 



a J 2tt 



(23) 



Using Eq. (j2"3"|) and Eq. (j2"2"j) and Carrying out the inte- 
gration first over the space and then over u>, one obtains 




According to Eq.© and Eq.© we have for the NLO 
correction 



(i) 



E, 



where Ejp 



lim 

b/a — >o 



lim 

L/b^a 



{Ei 1 



E; 

E, 



(1) 

6 

(1) 

bi 



2E, 



(i) 



(25) 



This computation is obviously complicated and plagued 
with a multitude of infinities. As explained before using 
the usual renormalization programs in conjunction with 
our box renormalization scheme, should eliminate all of 
the infinities, as might be apparent from the above equa- 
tion. However, proper regularization schemes should still 
be implemented and proper care taken when handling 
these infinite expressions. For example, the summation 
appearing in the squared form in the first term of the 
last part of Eq. ([24]) is infinite. We want to use the Abel- 
Plana formula to convert this sum into an integral. How- 
ever this sum does not satisfy the stringent requirements 
stated in the Abel-Plana theorem for such a conversion 
4ol 41 1 . However our box renormalization scheme pro- 
vides a solution: We first expand the square as a double 
sum. Then we subtract these double sums as indicated in 
Eq. (|25p . Now we can expect this new summand to sat- 
isfy the requirements for the Abel-Plana theorem. Then 
all the infinities actually cancel and the result for the two- 
loop correction reduces to (see Appendix for details): 



-Att 2 



/(a) -/(&) + 2/( 



L — a, ,L — b , 2 /„, , , . ,L - a, ,L — &, 

— -) - 2/(— - ) + - B(a) - B{b) + 2B(^—) - 2B(——) 



ds 



(26) 



r 



where, 



/(a) = B(a) 



B(a) 



I 

am 



coth(am) 
Am 



and B(a), defined by the following expression 
B(a) 



1 



t 2 7T 2 



m* 



(27) 



(28) 



refers to the so called branch-cut term in the Abel-Plana 
summation formula and is a finite quantity. Note that 
the last integral in Eq. (|26|) seems to diverge so it must 
be properly regularized, and this crucially depends on our 
box renormalization program, as we shall explain below. 



However, before we engage in this calculation, we want 
to raise an important point: If we were to use the free 
counterterm in Eq. (|21[) . as is routinely done, this term 
would be absent, in addition to some minor differences. 
Therefore, one would easily obtain finite results which we 
like to dispute. Now let us proceed with our calculations. 
We prefer to use a regularization scheme for this integral 
term which is analogous to the zeta function regulariza- 
tion for the sums. That is, in that expression we set the 
power of the integrand to — i + a for the first two terms 
and — i + a' for the remaining terms. In the final stage 
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we let a and a' approach zero. Hence we will have 
B(a)-B(b) 



'0 

Hto)-B(M) 



(1 



B(a) - 5(6) 



r(-a)+2 



g(jg)-B(jg) 



For a and a' sufficiently small, this expression becomes 



2 B(a)-B( b ) (_± +ln2 



+4 



B(^)-B( 



In 2 ) (29) then integration yields 



-1 

the infinities cancel 



Now, if — = -2 

' a S(a) - B(b) 

The cancelation of these divergent quantities without any 
residual finite terms is the usual prescription in regulated 
theories and this is what we have used . Therefore the 
term in question becomes, 

| In 2 (B(a) - B(b) + 2B(^) - 2B(^)) . (30) 

This result is obviously finite, and we believe it could not 
have been obtained with any regularization or analytic 
continuation schemes in common use, other than our box 
renormalization program. Thus, Eq. (|26p becomes 



-Avr 2 



, L — a 



e/, L - ^ 



f(a)-f'(b)+2f'(-— )-2/'(-— ) 



, (31) 



1 One may argue that ambiguities always exist in problems where 
one has to subtract infinite quantities, and the Casimir problems 
certainly fall into this category. Two methods are in common 
use: First is the analytic continuation techniques which, although 
usually yield correct results, do not have a very solid physical 
justification and also sometimes yield infinite results. Second is 
the regularization schemes, which is what we have used. In the 
latter category when the problem is regularized, one can make 
a systematic expansion of the quantities in question in terms 
of the regulators. Then the terms which tend to infinity when 
the regulators are removed and the finite terms naturally appear 
separately. See for example Eq. (129 D . What is almost invariably 
done is to adjust the regulators so that the singular terms ex- 
actly cancel each other, i.e. without extracting any extra finite 
piece from the difference between the infinite quantities (see for 
example [3). This is also apparent in the leading term for the 
Casimir energy in Eq. ([9} where, as explained in the Appendix, 
The four changes of variables are equivalent to choosing four dif- 
ferent cutoffs. One could have adjusted them so that as usual 
the infinities cancel, but any finite term would remain. How- 
ever, the well known answer is obtained only when there is no 
remaining extra finite term in this subtraction scheme. This is 
the prescription that we have used. However, we do believe that 
this is a subject that needs further study. 



where f'(a) = /(a) + 



21n2 



B(a). This is the two-loop 



radiative correction for the work done on the plates (or 
two points in this case) while moving them from (-^, §) 
to (-if, §)■ Now, in order to compute the Casimir energy, 
proper limits must be taken, as indicated in Eq. (|25[> . 

Two particular limits are interesting to calculate. First 
is the large mass limit. To calculate this limit it is conve- 
nient to make the following expansion in the expression 
for B(a), Eq. 



1 



2a 



B(a) - — V^o(2amj) 

7T ' 



am>l 



-2am j 



2 y / 2~amj 

Using Eq. ([33]) and Eq. ([31]), Eq. (25]) gives 
-A In 2 



(32) 



(33) 



-2bn 




In the small mass limit it is easier to rewrite an ex- 
pression for B(a), such that its integrand appears in di- 
mensionless form, 

2a r°° 1 dt 
n 



B(a) 



(35) 



», i \ r 

Then by expanding the integrand in Eq. (|35[) one finds 

_L_ a f°° dt 1 

2m 7r 



B(a) 



1 



2m 



— aS. 



i W 

Note the explicit appearance of infrared divergences in 
this equation which is a generic feature of these problems 
in 1+1 dimensions [35}. In this limit the first term in 
Eq. (|31|) due to region al, for example, becomes 



/ 1 



\2am 
= a(S 2 



1 



2am 



1 2 
S ) + - In 2 

am 7r 



1 



21n2 



S) 



In 2 

■Km 



4(am) 2 
(36) 



Hence, taking into account the analogous contributions 
from the other regions, i.e. bl, a2 and b2, Eq. ([2"S"]) gives, 



as 



0. 



(37) 



This shows that NLO for the Casimir energy in the mass- 
less limit is zero. Most importantly the infrared diver- 
gences have also cancel completely using our regulariza- 
tion program. This is in sharp contrast to the analogous 
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result that can be extracted from Refs. 2^, 25| which is Eq. (|23|) . Now in the Eq. (j22f . after space integration, 
infinite. carrying out the summation yields, 



C. The Massless Case 

In the massless case it is sufficient to set the pole of 
the propagator to zero, i.e. one can set m = in the 

I 



EV>(c 



—A 
8a 



du> 



au> coth(acj) 



duj'duj 



lo — uj 2 — au'uj coth(aa/) + au) u> coth(aa;) 
2co' 2 lu 2 (lu' 2 -co 2 ) 



(38) 



where, 



Pi = 



dp 



pcoth(p) — 1 



. , . p' 2 - p 2 - p'p 2 coth(p') + p' 2 pcoth(p) 
P 2 = / dp dp 



a 12 2/ 12 2\ 

4p P (p -P ) 
p = au and p' = aw'. 



Therefore, 



Ei 1 



n L — a L — b 
a -6 + 2 2 



(39) 

P2 = 0, (40) 



since Pi and Pi are independent of a, 6 and L. This 
result is in exact agreement with the small mass limit 
calculated in previous subsection. In figure we show 
our results for the zero and first order Casimir energies 
for the massive and massless cases. 

Note that we are explicitly assuming S m —o(x) ^ 0. Al- 
though this does not happen to make any difference in 
1+1 dimensions, using our prescription, we like to stress 
that this quantity should not in general be a priori set to 
zero. This is in contrast to the view expressed in for ex- 



ample Refs.[3j||33,|3a|39|]. This is yet another important 
counterexample for the validity of analytic continuation: 
As is well known the massless limit of the analytic con- 
tinuation of the mass counterterm in (/> 4 theory is zero 
for space-time dimension bigger than two. However one 
cannot renormalize the massless theory without the mass 



counterterm (see for example [34j). 



IV. CONCLUSIONS 

We have introduced a new concept in this paper. We 
have insisted that the renormalization program should 



-0.01 
-0.02 
-0.03 
-0 .04 
-0 . 05 



-0.06 




^10=0 



/ 



/ 



FIG. 2: The Casimir energy and its first order radiative cor- 
rection for the m = 1 and m = cases as a function of 
distance between the plates a. Note that E^L Q = 0. We have 
used the following conventions H = 1, c = 1, and A = 0.1m 2 . 



completely and self-consistently take into account the 
boundary conditions or any possible nontrivial back- 
ground which break the translational invariance of the 
system. We have shown that the problem is self- 
contained and the above program is accomplishable. To 
be more specific, there should be no need to import coun- 
terterms from the free theory, or even supplementing 
them with the attachment of extra surface countertcrms, 
to remedy the divergences inherent in this theory. In 
general this breaking of the translational invariance, is 
reflected in the nontrivial position dependence of all the 
n-point functions. As we have shown this has profound 
consequences. For example in the case of renormalized 
perturbation theory, the counterterms and hence the ra- 
diative corrections to parameters of the theory, i.e. m and 
A, automatically turn out to be position dependent. In 
this regard we disagree with authors who use the former 
counterterms (see the Introduction for actual references). 
Obviously we still need a few experimental input param- 
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eters for the complete renormalization program, such as 
m(x) for some x. However, the interesting point is that 
the theory then completely determines m(x). 

Secondly we have used a supplementary renormaliza- 
tion scheme, which is originally due to T.H. Boyer 3l|, 
along side the usual renormalization program. In com- 
putations of these sorts, there usually appears infinities 
which can sometimes be removed by the usual renormal- 
ization programs that often contain some sort of analytic 
continuation. These procedures are sometimes ambigu- 
ous. Our scheme is simply to confine the whole physical 
system in a box, and to compute the difference between 
the values of the physical quantity in question in two dif- 
ferent configurations. Use of this procedure removes all 
of the ambiguities associated with the appearance of the 
infinities, and we use the usual prescription for remov- 
ing the infinities in the regulated theory. Moreover all 
of the infrared divergences cancel each other. Using our 
method, we have computed the zero and first order radia- 
tive correction to the Casimir energy for the massive and 
massless real scalar field in 1+1 dimensions. For the zero 
order, our results are identical with what exists in the 
literature. However, our first order results are markedly 
different from those reported in Refs.0, 25 1. Our re- 
sults for the massive case is different from theirs due to 
the aforementioned conceptual differences. Moreover, we 
disagree with their results for the massless case obtained 
as the limit of the massive case, which is infinite in 1+1 
dimensions. Our analogous result is zero. The authors 
refer to the "exact" results obtained in Refs. 3y, l37[ as a 
verification of their massless limits in 3+1. However, in 



the latter references the authors effectively set 5 m equal 
to zero in their massless cases, or its equivalent. This 
is our second main difference in approach to the prob- 
lem. As mentioned before, we believe that 5 m should 
not be arbitrarily set to zero even in the massless case, 
since in that case the renormalization conditions can no 
longer be fully implemented, although the theory is still 
in principle renormalizable. 
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APPENDIX 

In this appendix we present the details of the calcu- 
lations leading to Eq. (1251). T he Abel-Plana summation 
formula (see for example [41|) is: 



/(it) - f(-it) 



2irt 



(It. 



E/(n) = -4r+f° mdx+i f 

n=l Z J ° Ja 

(A.l) 

In order to obtain first-order radiative correction, just 
as we discussed in section [2] for leading term, we need 
to compute — Ef \ Using Eq. (|2"4"|) . we get 



J 



E^ - E?> = 



(i) 



-Att 2 



i coth am — 1 bm coth bm — 1 (L — a) coth 



(L~a)rn 



E 



Aam 2 



46m 2 



2(L - a)m 2 



-S(a,n)S(a,n') - rS(b, n) S(b, n') + - 4 S( — - ,n)S( L - ,n' 
a b L — a 2 2 



1 (L-b) coth 



(L-b)r 



2(1 - b)m 2 
4 L-b . Qf L-b , 



(A.2) 



where S(a, n) = ( m 



2^2 \ -i/2 



Using the Abel-Plana formula Eq. (|A.1[) , and simple changes of variables in the 
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integrals, we obtain 



^(-q — A7r 2 I am coth am — 1 6m coth 6m — 1 (L — a) coth ^ L 2 °^ m — 1 (L — 6) coth 



(L-b)m 



4am 2 



46m 2 



2(L - a)m 2 



2(L- 6)m 2 



S(n,n) | S(6,ra) 2S{^,n) | 2S(^,n) 
2am 26m (L — a)m (L — b)m 



L — a 



, L — a 



+ - 5(a, n) - 5(6, n) + 2S(— - , n) - 2S(— - , n) 



f/.s' 



o %/ m 2 + s' ; 

' L-b\Q( L-b 



B(a)S{a,n) 5(6)5(6, n) 4B(^2)5(^2, n) 4B(^)5(^,n) 



L — a 



L-b 



Using Eq. (|A.1[) again, and making appropriate changes of variables to make the integrals dimensionless, all the actual 
infinities cancel and we finally obtain 



-Att 2 



f(a) f(b) + 2/(^) - 2/(V^) + I (B{a) - B(b) + 2B(^) 2b£-±) 



(A.3) 



It is important to note that all these cancelations are 
easily accomplished using our supplementary box renor- 
malization scheme. On a minor note, it is interesting 
to note that the changes of the variables leading to the 
cancelation of infinities are, surprisingly, equivalent to 



setting different cutoff rcgularizations on the upper lim- 
its of the integrals. Equation (|A.3j) is our main equation 
for the NLO Casimir energy, and appears in the text as 
Eq. ([26]) . and is analyzed further there. 
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